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Abstract 

Debye screening of static chromoelectric fields at high temperature is in- 
vestigated at next-to-leading order through one-loop resummed perturbation 
theory. At this order the gluon propagator appears to give rise to strong de- 
viations from a Yukawa form of screening. Generally, an oscillatory behavior 
is found which asymptotically becomes repulsive, but in a gauge-dependent 
manner. However, these features are strongly sensitive to the existence of 
screening of static magnetic fields. It is shown that a small magnetic screen- 
ing mass can restore exponential screening with a gauge independent value 
of the screening mass, which depends logarithmically on the magnitude of 
the magnetic mass. Recent results obtained in temporal axial gauge, which 
instead indicate an asymptotic (repulsive) power-law behaviour of screening, 
are also critically discussed. In order to arrive at a gauge-invariant treatment 
of chromoelectric screening, Polyakov loop correlations are considered, both 
with and without dynamical gauge symmetry breaking. Again a crucial sen- 
sitivity to the scale of magnetic screening is found. A detailed comparison 
of the perturbative results with recent high-precision lattice simulations of 
the SU(2) Polyakov loop correlator is made, which are found to agree well 
with the perturbative result in the symmetric phase when a magnetic mass 
~ g 2 T/A is included. 



I. INTRODUCTION AND SUMMARY 



It is well established that quantum chromodynamics (QCD) at sufficiently high tempera- 
ture and/or density is in a deconfined phase. Although the coupling g remains uncomfortably 
large up to astronomically large energy scales, it is hoped that this regime may be accessible 
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by perturbative quantum field theory at finite temperature and density. Indeed, perturba- 
tion theory works reasonably well at the energy scale set by the temperature, but it runs 
into infrared singularities when probing the softer scale gT, where a quasiparticle picture 
becomes relevant. 

The leading-order dispersion laws of the quasi-particle excitations are determined by the 
high-temperature limit of one-loop Feynman diagrams, the so-called "hard thermal loops" 
(HTL) which can be understood in classical terms But already at next-to- 

leading order the dispersion laws receive contributions from all orders of the conventional 
perturbation series. It has been shown in particular by Braaten and Pisarski [[J that an 
accurate and gauge-independent calculation of corrections to the HTL dispersion laws []5],|6| 
requires an improved perturbation theory which resums HTL contributions. This need is 
generic in perturbative thermal field theory [0,0, but in nonabelian gauge theories there is 
another, pernicious barrier for perturbation theory. Static magnetic fields are not screened at 
the HTL level, and the self-interactions of these lead to a breakdown of perturbation theory 
at a certain loop order depending on the quantity under consideration. The corresponding 
infra-red singularities are commonly expected to be cured by the dynamical generation of 
a magnetic screening mass ~ g 2 T, but its nature is still unclear. By superficial infra-red 
power counting, the next-to-leading order corrections to the dispersion laws are still below 
this critical order, however mass-shell singularities tend to generate a logarithmic sensitivity 
to the magnetic mass scale . 

In particular, such a sensitivity to the magnetic scale is found in a next-to-leading or- 
der calculation of the chromoelectric screening mass ||11|| , which will be recapitulated in 
Sect. 2. Such a nonabelian analogue of the classical Debye mass is supposed to account for 
an exponential screening of colour charges, and this will presumably provide an important 
characteristics of the hypothetical quark-gluon plasma |T2[ . There may however be substan- 
tial corrections to the pre-exponential part of the screening function, or even non-exponential 
behaviour on sufficiently large distances. For instance, in an electron gas it has been found 
T3fl that quantum corrections give rise to an asymptotic 1/r 10 behaviour at order h 4 , which 



eventually supersedes exponential screening. 

The main objective of this paper will be the discussion of the various possibilities to 
define screening functions and of the results found at next-to-leading order. The simplest 
possibility is to inspect the chromoeletric field induced by a weak static colour source, 
which is however a gauge dependent quantity. In Sect. 3 the covariant gauge result is 
analysed and, if taken at face value, it indicates a strong departure from the expected 
Yukawa-type potential. The pre-exponential screening function oscillates and approaches a 
negative value asymptotically, signalling a repulsive exponential tail. However, the analytic 
structures which gives rise to these phenomena turn out to be strongly sensitive to the 
existence of a magnetic screening mass, and are in fact found to be largely tamed by the 
latter. A small magnetic mass can restore exponential screening with a gauge independent 
Debye mass. In Sect. 4 recent results [0,0 obtained for the temporal axial gauge are 
discussed, which qualitatively differ from the covariant gauge results and which seem to 
imply an asymptotic repulsive power-law behaviour in place of exponential screening. The 
motivation for using the notoriously troublesome temporal gauge at finite temperature is 
its relation to the correlation of two chromoelectric field strength operators. Evaluating 
the latter in covariant gauge at next-to-leading order reveals a gauge dependence which 
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makes it clear that one cannot attribute direct physical meaning to it either. Moreover it is 
argued that the different analytic structure which is responsible for the different asymptotic 
behaviour found in Refs. ||14| , |T5|| may depend on the prescription for the additional poles of 
the propagator in temporal axial gauge. In Sect. 5, the manifestly gauge invariant Polyakov 
loop correlation (PLC) is used as a basis to determine the screening function. As has 



been most recently shown in Ref. |L6|, the large-distance behaviour of the PLC is strongly 



sensitive to the existence of a magnetic mass, and its introduction leads to exactly the same 
value of the Debye screening mass as obtained from the gauge-independent pole of the gluon 
propagator. In Sect. 6, next-to-leading order screening is discussed for the scenario of a 
static A condensate. While perturbation theory involves even more uncertainties in this 
case, it indicates a stronger sensitivity of the electrostatic potential to the magnetic scale. 
In Sect. 7 recent lattice results on nonabelian Debye screening are reviewed which seem to 
favour enhanced exponential screening at moderate distances. The data are found to agree 
well with the results for the PLC in the symmetric phase when a magnetic mass m m g 2 T/4 
is included. 



II. A PROPAGATOR-BASED DEFINITION OF THE DEBYE MASS 

In linear response theory, the chromoelectric field (E l (x)) induced by a single external 
source J is fully determined by the gluon propagator, without the need to consider higher 
vertex functions. With only one source there is also only one direction in colour space 
so that even the nonabelian field strength operator is linear in the gauge potentials (the 
commutator terms vanish trivially), and also covariant conservation of the external current, 
(D^A)).] 11 = reduces to simple transversality, <9 M J M = 0. In particular, the longitudinal 
electric field in momentum space is given by |Ti 



where 



h l K 2 



n L = n 00 - — iiy - -^-iv, (2) 

with Tlav the gluon self-energy which is diagonal in colour space in the absence of gauge 



symmetry breaking. (We follow the conventions of Ref. JT7| and denote 4-vectors by capital 
letters.) 

In the static case, Eq. (|j) reduces to ITl = n o and the electric field induced by a static 
source J M (x) = Q5q6 3 (x) is determined by 



<^(x)> = ^(W) (3) 



with the potential $ 



$(r) =Qj 



d 3 k 



e 



kr 



(27r) 3 P + n 00 (£;o = 0,A;) 
Q r°° e ikr - e~ ikr k dk 



(2ti) 2 J-oo 2ir fc 2 + n 00 (0, \k\)' 



(4) 
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If Iloo(0, k) is an even function in k, the last two pieces in Eq. (|4]) can be evaluated by 
closing the contour in the upper and lower half plane, respectively. At leading order in the 
high-temperature expansion, II 00 (0, k) is just a constant, 



n 00 (o, k) 



e 2 T 2 



2 T 2 



rrir 



with e 2 = (N + Nf/2)g 2 for colour group SU(A^) and Nf flavors, and Eq. 
simple poles at k — ±im , yielding 



Q 

— < 

Airr 



-mor 



(5) 

involves just 
(6) 



The next-to-leading order contribution to is down by one power of g rather than 



g 2 because of the "plasmon effect" [pl|. Whereas in Abelian theories Un U is a manifestly 



gauge independent object, in the nonabelian case it will generally depend on the gauge fixing 
parameters. Indeed, the infrared limit of n o turns out to be gauge dependent at the order 
gml, viz. 



19 



n 00 (0,A;^0) 



1 + a — 9 , 



where a is the gauge parameter of covariant gauges and 

N> 



1 + N f /(2N) 



(7) 



(8) 



However, this general gauge dependence does not mean that the screening function $(r) 
as defined in Eq. (|J) is altogether unphysical. If {k 2 + n o(0, k))" 1 still has simple poles at 
k = ±im one can prove on an algebraic level that their position, and therefore m, is gauge 
fixing independent [fj(]]. Hence, the screening mass m in the exponent of Eq. (|6|) can be a 
physical quantity, while the pre-exponential factor will depend on the gauge choice. This 



leads to the self-consistent determination of the Debye mass through |11 



m 



n 00 (o, k) 



k 2 = 



(9) 



rather than n o(0, k — > 0) which is usually taken as its definition ||18|| . 

The identification of n o(0, k — > 0) with the screening mass is in fact deficient already 
in the Abelian case. The infrared limit of n o is directly related to the second derivative of 
the thermodynamic potential with respect to the chemical potential [^TJ , and from this one 
knows |TI 



2 T 2 / o 2 

noo(o,fc-o) = l=-[i-^ + 



V3e 3 

47T 3 



+ 



(10) 



However, this result is not renormalization-group invariant. This is repaired by adopting 
(0), which amounts to adding 



n 00 (o,AO| 



-n 00 (0,A:^0) 



e 2 T 2 / 2& 



9tt 2 6tt 2 



In — + 7i; 



+ 0(e 4 
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where p is the mass scale introduced by dimensional regularization in which minimal sub- 
traction has been performed. The coefficient of the logarithmic term in ( |TT1) is exactly such 
that dm 2 /dp = because de/d(lnp) = (3(e) = e 3 /(127r 2 ) + 0(e 5 ). 

Turning again to the nonabelian case, it is now clear that we need more than only the 
infrared limit (|7|) of the next-to-leading order gluon self-energy. Since we are only considering 
the static case here, it is in fact possible to give the complete next-to-leading order result 
for rioo- In general such a calculation would involve the rather complicated propagators and 
vertices of the Braaten-Pisarski resummation program. However, as shown in Refs. [22]|§, 
for just the next-to-leading order contribution in static Green's functions, this resummation 
scheme boils down to the simpler ring resummation of Gell-Mann and Brueckner |23| . There 
one has to keep only the zero-mode contributions in the sum over Matsubara frequencies, 
and resummation consists only of the inclusion of the Debye mass in the longitudinal gluon 
propagators. 

The static ring-resummed propagator in general covariant as well as Coulomb gauge 
(with gauge parameter a) reads 



A 



fJ,V 



p =0 



p2 _|_ m 2 M 



a 



P P 



(12) 



JPo=0 



P \ P ) IP 

In these gauges, the complete next-to- leading order contribution to n o(0, k) is found as |TT 

d 3 ~ 2e 



5Tl o(k = 0,k) = gm 
2(m 2 - k 2 ^ 



p 



(2vr; 



3-2e 



p2 _|_ m 2 



+ 



p (q + m? 



+ (a - l)(k 2 + m 



2) p 2 + 2pk 
p 4 (q 2 + m? 



(13) 



where q = p + k. (Here dimensional regularization has been used when separating the 
static modes from the sum over Matsubara frequencies |22 ; the limit e — > gives a regular 
expression because of the odd integration dimension.) 



The new definition (||) requires to evaluate at k 2 



-m 



There the gauge dependent 



piece proportional to a vanishes algebraically, before doing the integrations, but the integral 
itself is linearly singular on the "mass-shell" k 2 = —m 2 . However, introducing a small 
infrared cutoff and taking the limit k 2 — > —in 2 before lifting the cutoff removes this term 
completely (cp. p3] , p5|j ). This can be achieved in a gauge invariant way either by dimensional 
regularization or by taking the symmetric limit of the Higgs mechanism |26| . 

The third term of the integrand in (^) is logarithmically singular as k 2 — > — m 2 . This 
singularity is caused exclusively by the massless denominator in the spatially transverse part 
of the gluon propagator ([T2]). A magnetic screening mass m m would screen this singularity, 
and because the latter is only logarithmic, the coefficient of the corresponding logarithm is 
unambiguously determined by (p~3|), 



5n 00 (o,fc) 



k 2 - 



g 2 NmT m 
In 



2tt 



(14) 



up to terms that are regular as m m — > 0. Assuming that m m ~ gm, the next-to- leading 
order contribution to m 2 is found to be of order ghi(l/g) rather than g, 
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Sm 2 s/W , 1 . 

7^ = — 9l % + 0(9) ' (15) 



which is positive, at least at weak coupling g <C 1, contrary to expectation [JT8[ 



The sublogarithmic terms cannot be calculated completely, because the presumed phe- 
nomenon of magnetic screening is nonperturbative p7|p8| . However, in order to obtain 
an estimate of the full 0(g) contribution, let us assume that a simple replacement of 
1/k 2 — > l/(k 2 + m^J in the transverse part of the static propagator (]TJ) correctly sum- 
marizes the effects at k ~ g 2 T. Then we may go on to evaluate the remaining contributions 
in (|T3|), which leads to 0, 



5m 2 v 7 ^/^ m / m 2 2m + m m 1 m m \ 2 

^^r = ^^^( (1 -4^ )ln ^ — 2-2^J +0( ^- (16) 

From the point of view of the original Braaten-Pisarski resummation programme || , the 
Debye mass m appearing on the right-hand side of the above equations as well as in Eq. (|9]) 
should be identified with the HTL value tuq in order to be consequently perturbative. How- 
ever, by the very introduction of a magnetic mass term we have already stepped beyond the 
resummation of HTL contributions. One could therefore equally well consider a resumma- 
tion of the correction terms to the classical value of the Debye mass at the dressed one-loop 
order, and determine m from Eq. ( jl6|) in a fully self-consistent manner. Diagrammatically, 
this corresponds to including not only "daisy" but also "super-daisy" diagrams |29| . 



Of course, all this does not affect the order 0(g) under consideration, but for realistic 
values of g the two possibilities will lead to noticeable differences numerically. The functional 
form of Eq. (|16j) is such that a self-consistent evaluation will generally give larger corrections 
5. We shall later adopt the variation caused by this as a measure of the uncertainty of 
the one- loop result, together with the more conspicuous uncertainty from the value of the 
hypothetical magnetic mass m m . 



III. CHROMOELECTRIC SCREENING FUNCTIONS IN CO VARIANT GAUGES 



Let us now ignore for the moment that the full function $(r) is gauge dependent and 
consider what screening behaviour it is purporting beyond the presupposed Yukawa form. 

Away from the singular points k = ±im, 5n o(0, k) as given by Eq. ( JJ~5| ) is regular so that 
there seems to be no need for the introduction of a magnetic mass as concerns the evaluation 
of Eq. (Q). Without a magnetic mass, n o(0, k) including the next-to-leading order is given 
by 



n 00 (o,fc) 



o gy/3N'm m 
m ° + 2~n 



m 2 — k k a — 2 

: — arctan 1 

mk m 2 



(17) 



If this expression is used to compute the full screening function $(r) from Eq. (Q), one 
finds a surprising behaviour [14]: $(r) decays exponentially, but the pre-exponential factor 
oscillates before approaching a negative constant for very large r. 

This strange result can be understood by inspecting the analytic structure of D^ 1 = 
(k 2 + Uoo(0,k)). In the left half of Fig. 1 this is displayed for the first quadrant of the 
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complex k plane (the others are given by reflection on the real and imaginary axes). At 
k/m = i, D^ 1 no longer has a simple zero, but a logarithmic branch singularity, and 
there is a branch cut from i to oo. The original zero of Dj}, however, still exists: it has 
moved to the right (and also left) of the imaginary axis. These complex poles of D L lead 
to an oscillatory behaviour of $(r); for the particular coupling constant chosen in Fig. 1, 
t = gv3N'/8 = 0.25, and gauge parameter a — 1, they contribute a term proportional to 
cos(0.313a;)e _1 ' 208:r /:r, where x = rm . There is, however, also the contribution from the cut, 
which adds a term —f(x)e~ x /x with a strictly positive function /, so that asymptotically, 
for very large x, the behaviour is that of a repulsive Yukawa potential with screening mass 
m . Both, the function / and the position of the complex pole, are gauge dependent. 




0.25 0.5 0.75 1 1.25 1.5 1.75 2 
Re k 



0.25 0.5 0.75 1 1.25 1.5 1.75 2 
Re k 



FIG. 1. The analytic structure of DJ^ = /c 2 +IToo(0, k) in units where 1 = uiq and with coupling 
t = g\/3N' /8 = 0.25. The left side corresponds to the case of zero magnetic mass, the right side 
to m in /m = 0.25. Full lines correspond to WleD^ 1 = 0, dashed ones to QmD^ 1 = 0, and their 
intersections to poles of Dl- The thick part of the imaginary axis marks the location of the branch 
cut; otherwise the imaginary part vanishes on both the real and imaginary axes. On the left, there 
is one (gauge dependent) complex pole; on the right, there are two poles on the imaginary axis: 
the lower one is at the gauge-independent position k = im, and there is a gauge-dependent one 
just below the branch singularity. 



However, if one allows for a small magnetic mass m m , the analytic structure, and therefore 
$(r), changes drastically. n 00 then becomes 



n 00 (o,fc) = m 2 + 



2 i gV3N'm 



2tt 
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— m m r 
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arctan • 



m^k ^ 
k 



m + m-n 



m m > . 



7 



where m m has been introduced in a gauge-invariant manner by mimicking the Higgs mech- 
anism [p6]. Choosing m m /m = t = 0.25, Eq. (^) gives a self-consistent Debye mass of 
m = 1.335m , and for these parameters the analytic structure of Dj} is rendered in the 
right half of Fig. 1. There is now a simple zero on the imaginary axis at k — im, and the 
logarithmic branch singularity has moved further up to k = i(m + m m ). The lines tfteDi = 
and ^stuDl = no longer intersect at complex values of k, but only on the imaginary axis. 
There is also a gauge dependent zero very close to the branch singularity, but the dominant 
singularity in is the one at k — im with m gauge independent. Thus $(r) no longer 
oscillates but decays exponentially for large x with gauge independent screening mass m. 

The strong dependence of the analytic structure of on the infrared behaviour of the 
transverse gluons of course means that perturbation theory cannot be trusted. Close to 
the branch singularity of , which arises at next-to-leading order, the correction terms 
become larger than the leading-order ones, and one clearly has to expect even more important 
contributions from higher orders. The appearance of gauge dependent poles is therefore just 
a manifestation of the incompleteness of the results in the vicinity of k = ±im Q . Assuming 
that all infrared singularities cure themselves in the complete results, e.g. by the generation 
of a magnetic mass, all gauge dependent singularities have to disappear according to the 
gauge dependence identities of Ref . [^U] . Whether in the final result the pole of the leading- 
order result survives at a shifted location cannot be established within the present resummed 
perturbation theory. However if it indeed does, Eq. (|l^) gives the leading order contribution 
to such a shift; the sublogarithmic terms of Eq. (|16|) on the other hand are just a simple- 
minded estimate. 

As we have seen, defining a chromoelectric screening function on the basis of the gauge- 
dependent gluon propagator allows one at most to extract a gauge independent exponential 
screening behaviour. All other details have to be considered unphysical. One longstanding 
proposal for doing better is to use a particular "physical" gauge fll8| , to wit, the temporal 
axial gauge. 



IV. NON-DEBYE SCREENING IN TEMPORAL AXIAL GAUGE? 



In recent work by Baier and Kalashnikov |T4| and, more rigorously following the Braaten 



Pisarski resummation program, by Peigne and Wong [IS], the screening function $(r) has 
been evaluated in temporal axial gauge (TAG), A = 0, and an asymptotic behaviour has 
been reported which differs qualitatively from the one obtained in covariant gauges. At very 
large distances $(r) ~ — 1/r 6 has been found, i.e. repulsive power-law behaviour, after an 
oscillatory regime at intermediate distances. 

The resummed calculation in TAG differs from the corresponding one in other gauges 
in that one cannot at once restrict to the zero modes of the resummed one-loop result. In 
TAG, the gluon propagator contains poles 1/Pq, which bring in contributions from resummed 
vertices, which would otherwise simply vanish in the static limit. These additional terms 
in n o(0, k) turn out to involve odd powers of k, i.e. a branch point at the origin of the k 2 
plane, and this is responsible for the asymptotic power-law behaviour. 

As mentioned in the introduction, a power-law asymptotic behaviour is known also from 
higher-order contributions to screening in a non-relativistic electron gas |L3l . More surprising 



here is the reversed sign, and also that such a behaviour should arise already at the next- 
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to-leading order, while the results in covariant and Coulomb gauges have a well-behaved 
Taylor series expansion in k 2 . 

Let us recall that the motivation for employing the rather awkward TAG is that there 
the electric field strength is linearly related to the vector potentials: E{ = —d^Ai. The 
longitudinal gluon propagator Dl, when computed in this gauge, is thus directly related to 
the correlation of two E operators. 

It has been argued [[H],[n|,[n| that the static correlation (El(x.)El(0)) is furthermore 



directly related to the free energy of a separated quark-antiquark pair in the singlet state and 
thus physically measurable. However, in this reasoning one has to consider external sources 
£ which couple to the electric field strength operator in order to construct a perturbing 
Hamiltonian if ext = J d 3 xi?"(x, £)£"(x, f). But E® is not a gauge invariant operator and 
the correlation of two E operators at separate points thus may depend on the gauge fixing 
parameters. Indeed, one finds that under a change of gauge condition f^A^ — > (/ M + Sf^A^ 
the correlation of two chromoelectric field operators varies according to 



5(E«(x)E e k (y)) = -gf abc J d 4 z(E){x)c c {x)c\z)5f,A d »{z)Et{y)) + (a,j,x~ e,k,y), (19) 

where c and c are Faddeev-Popov ghost fields, which make their appearance even in gauges 
which are otherwise ghost-free. 



d X % a d / \ a a ^lqjsT a 




f^*\ a 

+ X jyv*v\» +i/2 XX- wtv|A«M 

FIG. 2. One loop diagrams contributing to the correlation of two chromoelectric field strength 
operators. The wavy line with a dot is the ring-resummed gluon propagator and the dotted line is 
the Faddeev-Popov propagator. The crosses denote the extra vertices from the commutator term 
in the nonabelian field strength operator and the small box refers to the insertion of the Debye 
mass as a counterterm of the resummed perturbation theory. 

At the expense of having to keep the nonlinear terms in E, one can without further 
difficulty perform a resummed calculation of (El(— k)EL(k)) in gauges other than TAG. In 
order to obtain the next-to-leading order term, in almost any gauge except TAG one can 
restrict oneself to the zero-mode contributions, which more than repays the complications 
from the nonlinearity of E. A gauge which like TAG has a straightforward Hamiltonian for- 
mulation is the Coulomb gauge, which has the same zero-mode propagators as the covariant 
gauge. Evaluating the diagrams shown in Fig. 2 in Coulomb gauges with gauge parameter 
a, one obtains 



ab 



(El{-k)El{k)) = ^— 2 \1 - (20) 



9 



g NmT rm — k k a 

— arctan 1 H — 

mk m 2 



(27i)(k 2 + m 2 

g 2 NmT rm 2 - k 2 k 

H 77; — — arctan 1 

mk m 

k 2 + m 2 k " 

— arctan — 

mk m t 

g 2 NmT(k 2 + m 2 ) rk 2 - m 2 k 




87rk 4 



arctan hi 



mk m 
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The gauge parameter a does not drop out from this expression, which clearly shows that 
(El(x)El(0)) cannot be a physical quantity. It thus cannot be directly related to the free 
energy of two separated quarks, contrary to what is assumed in Ref. pO| , |3T|JT8| . 

In view of the non-Debye screening behaviour found in TAG, it is interesting to note that 
in covariant and Coulomb gauges {El(x.)El(0)) again has a series expansion in powers of k 2 
rather than k. This suggests that the different analytic behaviour found in TAG may well be 
an artefact of the treatment of the 1/Pq in the gluon propagator. In fact, such poles remain 
in the final expressions derived in Refs. |l4||n| , and are then evaluated by a principal- value 



prescription. 

But already at zero temperature, the principal- value prescription in axial gauges has been 
shown to be flawed in Wilson- loop calculations |33[, and at finite temperature the situation 
is even worse. Aq = is incompatible with strict periodicity in imaginary time [28], so 
when using the imaginary-time formalism (which is convenient for separating the zero-mode 
contributions) one has either to give up periodicity in the longitudinal propagator, which 
leads to rather complicated Feynman rules [j34}| , or one has to relax the condition A = 0. 
A regularisation of TAG through general axial gauges has been developed by Nachbagauer 
35f , and this procedure would again eliminate contributions from the resummed vertices, 



thus presumably removing the source of the contributions giving odd powers of k. And in 
the absence of odd powers of k, one would fall back to the task of finding the shift of the 
pole position, if any. If one defines an effective II ^ from the electric correlation (cf. Ref. JJI] ) 
as obtained in Coulomb or covariant gauges, Eq. (0), all the additional contributions are 
suppressed by powers of (k 2 + m 2 ), so that indeed the same Debye mass is obtained from 
Eq. (|9]) as through the conventional gluon self-energy. 

From all that it seems unlikely that a nonexponential asymptotic behaviour of Debye 
screening should occur already at relative order g. It might of course arise at higher orders, 
and in view of the largeness of g in actual QCD this would probably be an important effect 
already on not so large distances. 



V. GAUGE INVARIANT SCREENING FROM A POLYAKOV LOOP 

CORRELATION 

We have seen that the chromoelectric screening functions defined from the gluon prop- 
agator as well as from the correlation of electric field strength operators are strongly gauge 
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dependent in general and that only the exponential decay contributed by a singularity of the 
gluon propagator will be gauge independent and therefore of physical significance. If one 
is interested in the detailed screening behaviour beyond the value of the supposed electric 
screening mass, it is mandatory to first find a manifestly gauge invariant definition. 

A natural choice which can also be implemented without difficulty in lattice gauge theory 
is the Polyakov loop correlation (PLC) |36| . The Polyakov loop operator at the spatial point 
x is defined by 



1 / [P 



(21) 



where P denotes path ordering and r is the imaginary time. The correlation of two Polyakov 
loop operators is directly related to the free energy of a quark- ant iquark pair at the same 
separation |37|]. At lowest order (resummed tree-level), the connected part is given by 



(trfl + (r)trfl(0)) c - 1 



N 2 - 1 

8N 2 



' 'g 2 e- morS 
AurT 



+ 0(g 6 ) = + 0(g 



(22) 



and is given by diagram (a) of Fig. 3. 



1/2 OS. 




+1/2 & 




(a) (b) (c) 

FIG. 3. Resummed tree-level (a) and one-loop (b,c) diagrams contributing to the mean-square 
part of the Polyakov loop correlation. Here the wavy lines with a dot represent static longitudi- 
nal propagators while the transverse ones have been drawn as springs. Diagrams with vanishing 
next-to-leading order contribution have been omitted. Subtraction of the Debye mass countert- 
erm is understood to be included in diagram (b). With a finite magnetic mass, there is also a 
seagull- type self-energy insertion to be added to (b). 

The one-loop resummed correction to the mean-square correlation, fj£\ is given by the 
remaining diagrams of Fig. 3. Inspecting first large x = m Q r, one obtains 



/i 0) s +(1) //i 0) , - F(x) 



(23) 



1 + 



2tt 



1 7 
x( — In x — In 2 

v 2 2 



3. 
2' 



0(lnx] 



with 7 being Euler's constant and N' as defined in Eq. (§). The correction terms that are 
linear in x could be absorbed by a shift of the lowest order Debye mass ttIq — > m\ + 5m 2 , but 



there is also a term proportional to xlnx, which, if exponentiated, would lead to a potential 
that falls off faster than any Yukawa potential squared. Closer inspection of the Feynman 
diagrams 
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reveals, however, that the contributions oc x In x arise from singularities in 
momentum space which are strongly sensitive to the magnetic mass scale. 

Indeed, introducing again a small magnetic mass m m (see the appendix for details) 
changes the result (|23|) qualitatively for x ^> m/m m . The term oc xlnx contributed by 
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diagram 3b arises from a branch singularity coinciding with a double pole at k 2 = — m 2 . 
A magnetic mass m m separates these singularities and instead leads to a term linear in x. 
Its coefficient is given by the next-to-leading order correction to the Debye mass as derived 
from the gauge-independent pole of the gluon propagator in Sect. 2, 



X 



2vr 



— x 



(1- 



4m 2 



In 



2m + m n 



m T 
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2 



2m 



+ 0(x°) 



(24) 



The term oc x In x contributed by diagram 3c comes from a simple pole coinciding with 
a branch singularity at k 2 = —4m 2 . A small magnetic mass replaces the pole by two 
neighbouring branch singularities |l(J, yielding 



F (c) i 



x 



V3N>g 



2ty 



-(In — + 7) + 0(m/m m ) 



mi 



m 



(25) 



which grows only logarithmically at large x. 

In these results, the unphysical gauge modes have been made massive, too. If they were 
kept massless, they would contribute gauge-parameter dependent terms proportional to x 
to both F^ b ' and with opposite sign, leaving their sum unchanged. Absorbing all the 
linear terms of F(x) by a shift of the Debye mass implies F(x) —>■ F(x) + x5m 2 / (mm®) , 
thus leads to exactly the same next-to-leading order result for the Debye mass as from the 
gauge-independent pole of the gluon propagator, Eq. fli"6|) . 

In order to inspect the PLC screening function in more detail we shall need the form 
of F(x) also for moderate values of x. For x <C m/m m , the complete correction to the 
pre-exponential screening function F(x) of Eq. (^) reads [|3~6| 



F(x) 



1 H z x 



2tt 

-\a{2x) + 



■- lnx 
2 



ln2-2 + ^ 
2 2 



ln(2x) + 7 



+ 5 e 



l:r 



Ei(-4x) 



-e 2x m 



x 



-2x) 



(26) 



where Ei(x) is the exponential integral p8 ], x = mor, and 



dz 
z+l 



In 



z + 2 



(27) 



For large x, g(x) ~ ln(x)/x. 

Including the magnetic mass m m makes a difference for x ^ m/m m and we instead obtain 



= 1 + 



VW 1 
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27T m 
2 + /i 



(x) - -g(2x) + xh(2x) + |/i(4x) 



(28) 
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In 



2 + /x 



12 



where x = mr = xm/mo and \x = m m /m. We have further introduced 

h(x) = e x Ei(—x) 

and the function 




which we approximate by 

^^^riH^+i)- (so) 

2 + n z x 

In Fig. 4 the screening function F(x) exp(— 2xm/m ) is plotted logarithmically for the 
cases without and with a magnetic mass for the parameters used in Sect. 2. In order to get a 
notion of the uncertainties of the perturbative calculations, the next-to-leading order Debye 
mass m is evaluated both fully self-consistently (solving the nonlinear equation defined by 
Eq. fll6f) ) and strictly perturbatively (replacing m — > m in all quantities that are already 
of relative order g). One finds that for large x the dominant feature is exponential decay 
determined by the value of m. At small up to moderate distances the pre-exponential 
function is almost equally important. Notice that the effective Debye mass given by the 
slope of the tangent to the curves of Fig. 4 is generally smaller than the asymptotic value. 



log(S(x)) 




FIG. 4. The logarithm of the screening function Splc = F(x) exp(— 2xm/mo) of the Polyakov 
loop correlation for t = gV3N' /8 = 0.25, both without and with a magnetic mass. The 
short-dashed curve gives the case without a magnetic mass and m = thq, the full lines corre- 
spond to m m /mo = 0.25. The upper full line is for m from a fully self-consistent evaluation of 
Eq. (fi~6|), the lower one for a strictly perturbative one as discussed in the text. The long-dashed 
line marks the classical result, lnS* = —2x. 
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For the range of x covered by Fig. 4, the dependence of the screening function on the 
value of the magnetic mass is predominantly through the next-to-leading order value of the 
Debye mass rather than the explicit dependence of F(x) on m m . The latter dependence 
becomes relevant only for x ^> but there the perturbative result breaks down at some 
point, for F(x) also grows 1. 



VI. SCREENING WITH A CONDENSATE 

It has been speculated |39|.|4(| that the infrared divergences in the static sector of high- 



temperature nonabelian gauge theories cause the zero mode of Aq to develop a vacuum 
expectation value which breaks SU(iV) down to UtT)^" 1 . Nadkarni [40| has considered the 



effects of this scenario on chromoelectric screening in ring-resummed perturbation theory. 
From calculations in a unitary (diagonal) gauge, he suspected a strong modification of the 
Debye mass. 

In the diagonal gauge, where A only involves the diagonal generators of SU(N), correla- 
tion functions of the former correspond to the gauge invariant correlations of the eigenvalues 
of the Polyakov loop operator. The zero mode of A is restricted to v + e a T a where T a are the 
diagonal generators and v is the diagonal vacuum expectation value. In this parametrization, 
a nonvanishing v gives mass to the off-diagonal transverse static gluons A, which couple to 
the electrostatic fluctuations e through 



tr[A,v + e 



but there is no coupling between two electrostatic potential with only one transverse gluon. 
Since this was the only vertex in the next-to-leading order calculations in covariant and 
Coulomb gauges, it was maintained in Ref. [fHJ that all the results n o oc gm 2 were pure 
gauge artefacts. However, by inspecting the perturbative series in diagonal gauge, one finds 
that it ceases to exist for v — > 0. Denoting the mass acquired by the off-diagonal transverse 
static fields by ji = cg 2 T, where the value of c is not determined by perturbation theory, 
the loop expansion parameter turns out to be 1/c, which diverges in the symmetric limit. 
Hence, no statement can be made for the latter. 

Nevertheless, with c sufficiently large, the perturbation theory in diagonal gauge again 
makes sense and is in fact rather different from the case of no Aq condensate. The correlation 
of two e operators is manifestly gauge invariant by virtue of its relation to the Polyakov loop 
eigenvalues. In Ref. [[R]] this correlation was considered at resummed one-loop order only 
in the limit k 2 -C /i 2 . With the Feynman rules derived in Ref. j|(J it is straightforward to 
obtain the complete expression, which reads 

n 00 (0, k) = ml - £{l + 2 + - 2 + ^ I arctan - + -}. (31) 



For small k this coincides with the result given in Ref. 

The analytic structure of Eq. ( j3l|) is quite different from the ones we have encountered 
before. Because the correction term in Eq. ( PH ) is negative for large k and grows ~ fc 3 , 
there is generally a pole of ik 2 + n o(0, fc)) -1 on the real fc-axis. This is of course the 
dominant singularity for large distances and, taking the principal value, it would imply 
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an oscillatory behaviour ~ cos(Ar)/r. In fact, the residue is such that it even may start 
out repulsive at small distance, as shown in Fig. 5. However, this phenomenon is linked 
with the behaviour of the correction term in Eq. ( |31| ) at large values of k/fi, where it 
blows up due to the unfavourable large-momentum behaviour of the gluon propagator in 
the unitary diagonal gauge. Judging from Eq. ([H|), the new perturbation series breaks 
down for (k//i) 3 > c(m//z) 2 , which is in fact the region where the pole on the real axis 
arises. Since the other contributions to the screening function are suppressed at large k/fi, 
it seems reasonable to just exclude the contribution of the pole on the real axis while keeping 
the other analytic structures at smaller k. This can easily be done by a different choice of 
the integration contour in Eq. ([|). 



log(|S(x) |) 




FIG. 5. The logarithm of the screening function S(x) of the Polyakov loop eigenvalue correlation 
in the case of a nonzero Aq condensate with c = 1 and /i/mo = 0.25, where S(x) is normalized 
such that S c i = exp(— x). The short-dashed curve gives the oscillatory behaviour when the pole 
on the real axis is included (see text), which starts out with a reversed (repulsive) sign. The full 
curve is the result upon subtraction of the contribution from the real pole, and the long-dashed 
line marks the the classical screening function. 

With fi > mo/2 the only singularity is then the one at k « ±imo, so the asymptotic 
behaviour would be the classical one with only small corrections. However for fi < mo/2, 
which should be fulfilled at least with sufficiently small g, there is a branch cut on the 
imaginary axis superseding the original pole at k = ±imo. Hence, in this case, the magnetic 
mass scale is the dominant one in the asymptotic behaviour of the electric correlations. 
However, the imaginary part of (k 2 + Hoo(0, k))~ l along the cut is strongly peaked at k m 
±imo, so this dominance of the magnetic scale will occur only at very large distances. At 
intermediate distances the effective Debye mass is smaller than but close to the classical 
one, as shown in Fig. 5 by the plot where the real pole contribution has been subtracted. 

Altogether the resummed one-loop result for Debye screening with an A condensate 
appears more uncertain than with the conventional vacuum. The entire correction term 
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now is proportional to the coefficient c of the magnetic mass of the off-diagonal transverse 
gluonsQ. The value of c is inherently nonperturbative like the value of the magnetic mass 
which was invoked in the previous sections, but there the dependence was only logarithmic. 
A large value of c would provide a small loop expansion parameter of the new perturbation 
series obtained in the diagonal gauge, at least for sufficiently small momenta, but one can 
only expect c ~ 1. 



VII. COMPARISON WITH LATTICE SIMULATIONS 

Recently, high-precision lattice simulations of the Polyakov loop correlation in pure SU(2) 



gauge theory have been performed by Irback et al. f41fl . These simulations employ tempera- 
tures up to nearly 8 times the critical temperature so that one may hope that perturbation 
theory becomes applicable. Moreover, in SU(2) there exist also some (rather old) lattice 
results on the value of the hypothetical magnetic screening mass, so we can combine these 
results to test the various perturbative estimates for next-to-leading order correction to De- 
bye screening as presented in the previous sections, in which we had put in a magnetic 
screening mass by hand. 

There exist also some results for SU(3), both without ||2,33| and with quarks [33], but 
the corresponding lattice simulations have achieved less precision and involve smaller values 
of T/T c . Because of the latter, one has to expect more pronounced nonperturbative effects 
so that there is less chance for making contact with our perturbative considerations. Also, 
most of the lattice results for the magnetic mass, which is a crucial input in our perturbative 
estimates, have been obtained in pure SU(2). 



We shall therefore concentrate on the results given in Ref. |41| and compare them with 
our various perturbative estimates. The lattice simulations of magnetic screening in pure 
SU(2) gauge theory |^5[ have given m m = 0.27(3)g 2 T, which is consistent with a recent 



semiclassical result by Biro and Miiller [|4(| and also with a recent result from one-loop 



resummed gap equations in a non-linear sigma model for the magnetic mass term [27|. In 
what follows we shall plug in this value for the magnetic mass into our perturbative result for 
the next-to-leading order Debye mass ([16]) with the above error and evaluate it both strictly 
perturbatively and fully self-consistently. As discussed at the end of Sect. 2, in the strictly 
perturbative case we only resum the classical (HTL) Debye mass mo whereas in the second 
case we resum the full m, solving the resulting nonlinear equations for m numerically. The 
variation of the results will serve us as a crude measure for the uncertainties of the resummed 
one-loop calculation. 



In Ref. PTJ, the Polyakov loop correlation has been measured for lattice sizes 4 x 24 3 at 



g 2 = 1.43 and 4 x 32 4 at g\ = 1.33, at distances of up to 8 lattice units. In order to make 
contact with perturbation theory, the coupling constant has to be renormalized. This is done 
at short distances, that is, at the smallest available distance of one lattice spacing, yielding 
g 2 m 1.28 and g 2 m 1.19, respectively. The smaller value of the coupling, which pertains 



The magnetic mass of the diagonal gluons does not enter here and is in fact independent of the 
one provided by the assumed dynamical symmetry breaking. 
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to the larger lattice, corresponds to the rather high temperature T/T c « 7.8 according to 
recent work by Fingberg et al. |4q| . With this renormahsed value of g, Eq. ([16|) yields for 
m m /g 2 T = 0.27 ±0.03 

m_fl.22 T 0.02 , . 

m 11.36 =f0.04' 1 J 

where the upper and lower lines correspond to perturbative and self-consistent evaluation, 
respectively. (With the larger value of g the above values are almost unchanged.) 

A screening mass has been extracted from the lattice data by fitting the unsubtracted 
PLC to the mean-square behaviour of Eq. (]22|) at distances of 4 to 8 lattice spacings. For 
both lattices, this corresponds to values of x ~ 0.9 ... 1.8. Although these are not very large 
values, the mean-square approximation is probably quite good in SU(2), since there the 
next important correction term associated with 3 gluon exchange vanishes []3"E| ]. The results 
reported in Ref. [[11] are 



m 
m 




(33) 



i.e., a significant enhancement of the Debye screening mass over its classical value, which is 
even more pronounced for the smaller coupling corresponding to the higher temperature. [] 
However, for the moderate values of x ~ 1 ... 2 covered by the lattice calculations these 
numbers cannot yet be compared with the perturbative results which refer to the asymptotic 
behaviour. As we have seen in Sect. 5, it is important to take the pre-exponential screening 
functions into account when attempting to identify the asymptotic screening mass. From 
Fig. 4 it is clear that fitting the classical function (e~ mr /r) 2 on a restricted range of x will 
underestimate^ m. In Fig. 6 the lattice data are compared with the perturbative results 
of Sect. 5 for the screening function F(x) exp(— 2xm/mo). Fitting straight lines in this 
logarithmic plot gives an effective screening mass which reads 



m 
m 



■eff. _ J 1.16(9) PLC with m m and Debye mass from Eq. (32) ,„ 
0.97 PLC without magnetic mass and m = 



where the error gives the full variation caused by both the quoted error of the magnetic mass 
and the spread from a perturbative vs. self-consistent evaluation of the one-loop formulae. 
The PLC result without a magnetic mass thus leads to a slightly diminished effective Debye 
mass as compared with the classical value. On the other hand, the increase of the Debye 
mass according to Eq. (|i"6|) is not completely compensated by the curvature of F(x) and 
leads to an effective Debye mass significantly larger than the classical value, in accordance 



Debye screening in SU(2) gauge theory was also studied in Ref. [49], but at somewhat smaller 
distances and not directly in terms of the connected Polyakov loop correlation. 

3 A different approach to determining the Debye screening mass, which is based on the transfer 
matrix formalism, has recently been proposed in Ref. [50] and indeed leads to somewhat higher 
values than the conventional approach. 
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with the lattice data. As concerns the magnitude of the screening function, the lattice data 
lie in between the perturbative results with and without a magnetic mass. 




1.2 



1.6 



1.8 



FIG. 6. The logarithm of the screening function F(x) exp(— 2xm/mo) of the Polyakov loop 
correlation in the range of x = m^r covered by the lattice calculations of Ref . [^] . The long-dashed 
straight line gives the classical result. The short-dashed curve gives the one- loop resummed result 
for the mean-square contribution prior to the introduction of a magnetic mass. The full lines give 
the range of the results which include a magnetic mass m m = 0.27g 2 T. The upper and lower full 



line correspond to the upper and lower value of the Debye mass m in Eq. (32), respectively. The 
lattice data for the connected Polyakov loop correlation are given by the thick dots. The error bars 



as reported in Ref. [41| are negligible except for the last two data points. 



All of the above numbers depend crucially on the value of the renormalized coupling 
constant associated with the bare one on the lattice. So far we have followed the procedure 
of Ref. 
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where the former was determined by matching the short- distance behaviour to 
the tree-level result of the bare perturbation theory. Using instead the resummed tree-level 
result which includes the leading-order Debye mass necessarily gives a different renormalized 
coupling because one cannot go to smaller distances than one lattice spacing. In the case 
of the larger lattice and higher temperature, the latter procedure increases g$ = 1.33 to 
g\ ~ 1.53, while previously the renormalized coupling was diminished. With this larger 
value of gR, the results in Eq. (|34j) are only insignificantly increased, yielding now 



m 



eff- 



m 



1.18(11) PLC with m m and Debye mass according to Eq. (16) 
0.98 PLC without magnetic mass and m = m Q 



(35) 



whereas the lattice result ( |3"3"D is reduced to m/m = 1.00(4). Apparently this is in better 
agreement with the perturbative result prior to the introduction of a magnetic mass, but by 
directly comparing the complete screening functions with the lattice data, redone in Fig. 7, 
the case with a magnetic mass and an enhanced Debye mass from Eq. (|16"D , seems to be 
clearly favoured. 
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FIG. 7. Same as Fig. 6, but with a renormalized coupling of g 2 ~ 1.53, which is obtained by 
matching the short-distance behaviour to the resummed (rather than bare) tree-level result. 

Finally we briefly turn also to the case of an Aq condensate, where we have encountered 
a rather different analytic structure of the screening function. There the theoretical uncer- 
tainty is much greater, because the loop expansion parameter 1/c is always > 1 (with the 
above parameters 1/c ~ 4). In Sect. 6 we have found an oscillatory behaviour that started 
out repulsive, but we have argued that this result is not perturbatively stable. After discard- 
ing this part of the result, there remained a contribution from a branch cut at k = ±2im m 
instead of ±im. For very small coupling, this should be a distinctive feature, at least at 
sufficiently large distances. However, the lattice results involve only moderate values of x 
and furthermore, for the coupling constant that was used in the above lattice simulations 
and the assumed magnetic mass, we have 2m m /m = 0.7(2). For large x this suggests a 
diminished Debye screening mass, but at the not so large x for which lattice data exist, 
the effective screening mass may be different. Indeed, the screening function displayed in 
Fig. 5 has positive curvature (in contrast to those of Fig. 4), so for smaller x the effective 
screening mass now becomes larger. With the lattice parameters, the one-loop result yields 
m e ffJrriQ ~ 1.05, i.e. even a slight enhancement compared to the classical Debye mass, thus 
going in the right direction. However, as has been shown recently [^TJ, in the case of a non- 



vanishing A condensate the leading term in the PLC is not the mean-square contribution of 
Eq. ( P4) but one which involves only one power of the usual Yukawa-type potential. Fitting 
the latter behaviour to the lattice results of Ref. |^T[ gives m/m ps 3.4(1) (« 3.0(1) if the 
larger value of the renormalized coupling is used), which is far off from the perturbative 
result .0 



4 A single power of the Yukawa potential is in much better shape just above the deconfinement tran- 
sition [49 1, and it has been conjectured in Ref. [51| that there might be another, gauge-symmetry 
restoring phase transition not far above the deconfinement temperature. 



19 



In view of the largeness of the coupling and the uncertainty associated with the treatment 
of the magnetic mass, one certainly cannot expect next-to-leading order results to be really 
accurate, but the agreement of the lattice data with the screening function of the symmetric 
scenario with a magnetic mass m m ~ g 2 T/4 and the enhanced Debye mass according to 
Eq. (|IED is in fact surprisingly good. The range of x covered by the existing lattice cal- 
culations is however too small to really pin down the actual asymptotic behaviour. Even 
the qualitatively different asymptotics suggested by the TAG result, which shows oscillatory 



behaviour with a repulsive power- law tail [14], would manifest itself only for x considerably 



greater than 2. However, lattice simulations at larger values of x are probably prohibitively 
difficult, since the statistical errors increase very quickly with x. On the other hand, progress 
with (resummed) perturbation theory seems only possible through a better understanding 
of the nonperturbative magnetic sector of thermal nonabelian gauge theories. 
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APPENDIX: 

In this appendix some details to the evaluation of the PLC screening function are given 
when a magnetic screening mass is included, supplementing the calculations contained in 



the appendix of Ref. [36]. 



F^ b \ introduced in Sect. 5 and corresponding to the function /j of Ref. |]36||, is given by 

d 3 ke~ ikl 
(k 2 + m 2 ) 



= 8nre™ / ffi*"*' (*m a - OT 00 (0, kj) (Al) 



where <i 3 k = d 3 2e /c/(27r) 3 2e and 



5U 00 (k = 0, k) = g 2 NT I d 3 p ( - - - + — - 



m 2 m 



2(m 2 - ml/2 - k 2 ) , . ^ 2 , 2 . p 2 + 2pk 



+ (a - l)(k 2 + m 2 )-— -f— — - , (A2) 



with q = p + k. 
This leads to 



p 2 (q 2 + m 2 ) (p 2 + m 2 n ) 2 (q 2 + m 



Sm 2 g 2 NT((m + m m )r . . f , , (1 - e iVr ) 



z 2 + 2z + 2 - /i 2 /2 n-z 



+ - / dze~ mrz ^^ln 



vr Jc ' z 2 (2 + z) 2 2 + fI + Z 
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where the contour C encircles the positive real axis and \x = m m /m. 

The contour integral in (|A4[) receives contributions from a pole at z = and a branch 
cut starting at z = //, which thus is separated from the pole by the finite magnetic mass. 
The contribution from the pole reads 



Z7r[(l-^/4) 
and the one from the cut 



-mr In 



2 + fi n 2 + n 



+ ^ln^ 
4 2 + /x 



^(1-^/4) 
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(A4) 



(A5) 



Ei(-m m r) + 



A 4 Jmr 



Ei(-(2 + /i)mr)|. 



In general covariant gauge and with a magnetic mass, F^ c \ which corresponds to the 
function /// of Ref . [EfJ , reads 



with 



h 
h 
h 



F (b) = g 2 NT |_ 2 ^ _ Q j j2 + ( 47rr )2 e 2mr ^ _ ^ _ ^ j 
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(A7) 
(A8) 
(A9) 



The gauge-parameter dependent parts of exactly cancel those of F^ b '. The integral 
appearing therein yields 



cZ-J 



(l 2 + m 
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m/ 



J— (l - e-^) , 



(A10) 



so even separately these terms do not contribute to the linear terms in r that constitute 
the next-to-leading order shift of the Debye screening mass if the unphysical gauge degrees 
of freedom are made massive; if they are kept massless, I2 —>■ r/(8n), and these linear 
contributions cancel upon adding F^ b > and F^ c \ 

In Ref. |36| all of the above integrals were evaluated exactly in the case of a vanishing 
magnetic mass. With a finite magnetic mass, J 3 can still be calculated completely, which 
yields 
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Ii is unchanged by the introduction of the magnetic mass and reads 



[Aixry An 



with g as defined in Eq. (p7|) . 

J 5 , with a finite magnetic mass, can be brought into the form 



h = J rf 3 ke" ik - r (k 2 + 2m 2 - m 2 m /2)J(k) (A13) 



where J(k) is given by the double integral 



j(jfe) = J_ /' (A14) 
32tt^o ^^(l 

— x) + m 2 
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- 2x2/) 



[A; 2 + 4(m 2 + yH / )]i/A; 2 x?/(l — xy) + m 2 + yVF 



where W(k, x) = m^ + 2m m y Px(l — x) + m 2 . 
With m m = W = 0, 

]_ k 2 -\- 4m 2 

J(A:) = 167r 2 A; 2 (A ; 2 + 4m 2 ) ln 4m 2 ' (A15) 

which leads to 



e -2mr r r dze- 2mrZ (l +4Z + 2Z 2 ) 



I 5 = / ^_ _ i ln[-z(2 + *)] (A16) 

5 (4vrr) 2 8i7r 2 Jc z{l + z){2 + z) K ' 
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ln(mr) + 7 + 2e 2mr Ei(-2mr) + e 4mr Ei(-4mr) 



As has been pointed out first in Ref. [|16|1 , a finite magnetic mass removes the pole at 
k 2 = —4m 2 (corresponding to z = 0) and replaces it by neighbouring branch singularities at 
k 2 = —Am 2 and k 2 = —(2m + m m ) 2 . In the limit k 2 = —k 2 — ► 4m 2 , J(fc) yields 

1 2m — k 4m 2 — k 2 , , _. 

•J(*0 , , ~7 ^7 ln ^ ln — — , (A17) 

fc2^„ 4m 2 (167rmm m ) 2 4m 4m^ 

which leads to the large distance behaviour 

-(hi2m m r + 7 ), (A18) 



e 2mr m 



(47rr) 2 47rm 2 n 



up to terms that are down by powers of m m /m. 

Putting everything together leads to Eq. (p9|) , where the full function I 5 has been ap- 
proximated in terms of the function £(x) introduced in Eq. fl30|). 
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